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Inverse Problems for Systems of Variable Order Differential 
Equations with Singularities on Spatial Networks. 

V. Yurko 

Abstract. Variable order differential equations with non-integrable singularities are consid¬ 
ered on spatial networks. Properties of the spectrum are established, and the solution of the 
inverse spectral problem is obtained. 
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1. Consider a compact star-type graph T in with the set of vertices V = {uq, • • •, Vp} 
and the set of edges £ = {ci,..., Cp}, where Vi,... ,Vp are the boundary vertices, Uq is the 
internal vertex, and cj = [nj,no], ei fl ... n Cp = {uq} . Let Ij be the length of the edge ej . 
For each edge cj & £ we introduce the parameter xj G [0, Ij] such that Xj = 0 corresponds to 
the boundary vertices Vi,... ,Vp . A function Y on T may be represented as Y = , 

where the function yj{xj) is dehned on the edge Cj . 

Let rij , j = l,p, be positive integers such that rii > n 2 > ... > Up > 2. Consider the 
differential equations on T : 


y 


K), 


Xj) + 


1 - 1=0 




+ Qt^j(xj )) y)'"^(xj ) = Xyj(xj ), xj e {0,1 j), j = 1,p. 


Xa 


,(m) 


( 1 ) 


where A is the spectral parameter, qp_j{xj) are complex-valued integrable functions. We call 
^3 ~ ^l= o,nj -2 fhs potential on the edge Cj , and we call q = {qj}j^i^ the potential on 
the graph T. In this paper we study inverse spectral problems for system (1). We provide a 
procedure for constructing the solution of the inverse problem and prove its uniqueness. 

2 . Let ns construct special fundamental systems of solutions for higher order differential 
operators with regular singularities. Consider the differential equation 


n—2 

iy=yi^') + '^(^-^ + qj[x)^yO) = Xy, o: > 0 (2) 

3=0 ^ 

on the half-line. Let pi,..., /i„ be the roots of the characteristic polynomial 

n j-l 

^(h) = 'Yh = T^n-l = 0. 

j=0 k=0 


It is clear that pi + ... + /i„ = n{n — 2)/2. For dehniteness, we assume that pk — yj 7^ sn 
{s = 0, ±1, ±2 ,...) ; Re fix < ... < Re Hn , /ifc 7 ^ 0 , 1 , 2 ,..., n — 3 (the other cases require 
minor modihcations). Let 6n = n — 1 — Re{fin — hi) , <loj{x) = (lj{x) for a; > 1, and qoj{x) = 
qj{x)x'^^'^^^'^~^’^^ for a; < 1 and assume that qoj{x) G L(0, cxo), j = 0,n — 2. 

First of all, we consider the following differential equation without spectral parameter: 


n—2 

% := = y- 

3=0 ^ 

Let X = rexp(z<p), r > 0, p G (—7r,7r], x^ = exp(/i(lnr + iip)) and n_ be the x-plane with 
the cut along the semi-axis x < 0. Take numbers Cjo, 4 = 1,’^ from the condition 


JJcjo = (det[hj ^]j,u=T^) 
i=i 
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Then the functions 


Cj{x) = ^ Cjkx'"’", Cjk{x) = Cjo JJ A{fij + 


k=0 


are solutions of (3), and det[Cj‘^ ^\^)]j,u=T^= Moreover, the functions Cj{x) are analytic 
in n_ . Denote 


sn 


S=1 


-1 


(4) 


Sk = exp 


27ii{k — 1) 


n 


n I (z/ + l)7r 

, = \x ■. dxgx e [—,- 

\ n n 


Si — Sn-l, Sj. — Sn-2k+l U Sn-2k+2) k — 2,n', 


71 


Qk = \ x '■ axgx G max ( — tt, {—2k + 2) — ), min ( tt, (2n — 2A; + 2) — ) k k = l,n. 


n 


71 


n 


For X ^ SI equation (3) has the solutions ek{x), /c = l,n of the form 

(^k~^\^) = elexp{ekx)zku{x), z/ = 0,n- 1, 
where Zku{x) are solutions of the integral equations 

"I OO ^ ^ ^ 

Zky{x) = 1 + - exp{{ek - £j){t - x))^ Zkm{t)^ dt 

^ j=l m=0 

(here argt = argx, \t\ > |a;| ). Using the fundamental system of solutions {Cj{x)}j^i^ one 
can write 

n 

(ik{x) = ^ldlfj{x). (5) 

t=i 

In particular this gives the analytic continuation for ek{x) on n_. 

Lemma 1. The functions {ek{x)}x G n_ form a fundamental system of solutions of 
equation (3), and 

det[e^r^\x)],,^^^i^ = det[el-\^^i^. 

The asymptotics 

exp(efca;)(l + 0(a;“^)), \x\ ^ oo, x E Qk- (6) 

is valid. 

We observe that the asymptotics (6) holds in the sectors Qk which are wider that the sectors 
S^ . Next we obtain connections between the Stokes multipliers 
Lemma 2. The following relations hold 


Pkj = ^ij^kS j:k = l,n, 

n 

Wfdlj = ^ 0. 

i=i 

Indeed, for argx G (— tt,tt — —) we have, by virtue of (4)-(5), 


(7) 

( 8 ) 


ek{e"x) = '^ldlj{e^Y^Cj{x). 
i=i 


( 9 ) 
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It is easily seen from construction of the functions ek{x) that ei{e^x) = es+i{x). Substituting 
(5) in this equality and comparing the corresponding coefficients, we obtain (7). After that (8) 
becomes obvious. □ 

Now we consider the differential equation 


ioy = Xy = p'^y, a; > 0. (10) 

It is evident that if y{x) is a solution of (3), then y{px) satihes (10). Define Cj{x,\) by 

OO 

C,{x, A) = p-^^C,{px) = x^^ 

k=0 

The functions Cj{x,X) are entire in A , and det[C]‘^ ^\x,Xy\j^,^l^ = 1. From Lemmas 1 and 
2 we get the following theorem. 

Theorem 1. In each sector Sk^ = {p ■ argp G (^, ■^^^^^)} equation (10) has a funda¬ 
mental system of solutions Bq = {yk{x, p)}i.^Yn that yk{x,p) = yk{px), 


I") I 


X >1, z/ = 0, n — 1, 

(11) 

c \k,u=l~n 7^ 0) 

(12) 

, ^ 0, 

(13) 


\p\x' 


dei[y^^ ^\x, p)\k^^^TIl = D := det[i?^ ^]fc,^=T7i ^ 


i=i 

where the constant Mq depens only on Uj. 

The functions yk{x, p) are analogues of the Hankel functions for the Bessel equation. Denote 

C*{x, A) = det[C'^''^(a;, A)]^=o 7 t= 2 ;fc=M\n-i+n 

yj{x, p) = (^p^^-^Xn-2)/2^ det[y^^\x, p)]y=o^,k=Vti\j. 


Bkuipx) = 


R'(.exp{pRkx), \p\x>l, 


{px] 


fll-U 


\p\x < 1 , 


F^ipx) = 


exp{-pRkx), \p\x > 1, 

(pa;)"^“^“^", |p|a; < 1, 


UL{x,p) =yk\x,p){p’'Fk^{px)) \ U^’*{x,p) =yl{x,p){F^{px)) \ 


It ^11 

g{x,t,X) = '^{-l)^-^Cj{x,X)C*_j^^{t,X) = -^^yj{x,p)y*{t,p). 
j=i ^ i=i 

The function g{x, t, X) is the Green’s function of the Cauchy problem ioy—Xy = f{x), y^'^\0) = 
0, z/ = 0,n — 1. Using (11)-(13), we obtain 

\Uj)^{x,p)\ < Ml, \U^’*{x,p)\ < Ml, a;>0, p e Sk^, 


(14) 


dp 


dx‘ 


|Gf^(a;,A| < M^lx^^^-^l, 

n 

:g{x,t. A) < M 2 \px\ < Cq, t < x, 

i=i 


( 15 ) 


where Mi depens on Uj, and M 2 on Uj and Cq. 
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Now we are going to construct fundamental systems of solutions of equation (2). Denote 

n—2 

J{p) — ^ ^ Jm{,P)i 

m=0 


Jmip) — IpI 




HpI- 


t^" '^\qmit)\dt+\p\ 


m—n+1 


'\p\- 


\qm{t)\dt. 


Lemma 3. The following estimate holds 


^ n 2 

J(p) < IpI >1. Q j l9om 

I''! 


(f)| dt. 


Indeed, if Om — 'm < 0 , then Re{pn — Pi) > n — m — 1, and consequently 

p|p|-l POO POO 

JUp) < |pr-”+M / t^--nqm{t)\dt+ / \qUt)\dt) < \pr-+^ / \qoUt)\dt. 

^Jo J|p|-1 ^ Jo 

If 9m — m > 0, then 

i*00 POO 

Jmip) < |P|“' / \qmit)\ dt < |pr-"+' / |go™(t)| dt. 

do Jo 

Hence J(p) < (5|p|“^ |p| > 1, and Lemma 3 is proved. 

We now construct the functions Sjix, A), j = l,n from the system of integral equations 


Sf\x,\) = Cf\x,\)- j ■^gix,t,\)(^'^qmit)S^'^\t,\)'^dt, u = 0,n-l. 


n—2 


(16) 


m=0 


By (15), system (16) has a unique solution; moreover the functions Sj^\x,X) are entire in A for 
each a; > 0, the functions {S'j(a;, A)}j=xy;; form a fundamental system of solutions for equation 
(2), det[5'j''“^^(a;, A)}j.^=x;^ = 1, and 

S^J{x, A) = Oix'^^-’^), (Sjix, A) - Cj{x, \))x->^^ = o(a;^"-^i), a; ^ 0. (17) 

Let Sko^a = {p ■ P £ SkQ, IpI > a}, po = 2MiQ + 1. For k = l,n , p G SkQ^pQ we consider the 
system of integral equations 


^ 2 POO 

Ukuix, p) = U^jyix, p) + X] / Akumix, t, p)Ukmit, p) dt, X > 0, u = 0,n- 1, (18) 

m=0 do 


where 


Akvm ix, t, p) 


qmit)Fkmipt) 

pn-l-mFk^{px) 


k 

- p)F*ipt)Uf*it, p), 

i=i 

n 

E FApOUl{x,p)F'(pt)Up(t,p), 

j=k+l 


t < X, 

t > X. 


Using (14) and Lemma 3, we obtain 

poo Q 

P, / \Ak^,„(x,t,p)\dt < MiJ(p) <—^ 

I'’! 
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Consequently, system (18) with p G Sko,pQ has a unique solution, and uniformly in a; > 0, 

Uku{,x, p) - p) = 0(p"^), p e Sko,po- (19) 

Theorem 2. For a; > 0, p G SkQ,pQ there exists an fundamental system of solutions of 
equation (2), B = {Yk{x, p)}k=-i^ of the form 

^k\^^p) = p''Fku{px)Uky{x,p), 

where the functions Ukuix,p) are solution of (18), and (19) is true. 

The function Y^''\x,p) considered for each a; > 0, are analytic in p & SkQ,pQ, continuous 
in p & Sko,po ond det[y)!'‘^~^^(a;, p)];;. + 0(p“^)) as |p| —)■ cxd. The functions 

Yk{x,p) satisfy the equality 

Ykix^p) = ykix,p) - / (^^yj{x,p)y*{t,p)^(^^qmit)Y^"'\t,p)^ dt 

P j=l m=0 


n—2 


p' 


,n— 1 


p^) ( p)) 

j=k+l 


m=0 


Moreover, one has the representation 


Yk{,x,p) = ^hkj{p)Sj{x,\), ( 20 ) 

i=i 

where 

Kip) = ^°j(p)P^^(l + OipK). \p\ ^ oo, p G Sk,,p,. (21) 

The only part of the theorem that needs proof is the asymptotic formula (21). Let p be 
hxed, X < IpK- Then (13) and (20) become 


n \ 

USo(^. P) = Y ^). 

n 

Uko{x,p) = '^bkj{p){p)~^^x^^^->^Kj{x,X), 

j=i 

where 


( 22 ) 


Cj(x, A) = x-^tCjix, A), Sj{x, A) = x-'‘’Sj(x, A), Sj(0, A) = p(0, A) = c,o # 0. 

It follows from (22) that 

n 

Uko{x,p) - U^oix,p) = (bkjip)p~^" - bljp^^^->^^^x^^^->^Kj{x,X) 

i=i 

n 

+ -y(x,A)). (23) 

i=i 

Denote 

Fkiix,p) = Ukoix,p) - f/)?o(a;,p), 

Fk,s+iix^P) = - Fksitl, p)Ssix,X)cfQ^x^^‘’~^^‘^+\ s = l,n- 1. 
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Lemma 4. The following relations hold 

ibksip)p~'"^ - hl^p^^‘~^^)cso = J'fcs(a;,p), s = ifn, (25) 

s-l 

Tt.{x,p) = {(Ua{x,p) - UUx,p)) - Y^(h,j{p)p-^‘'-h%p^>-^‘')x‘‘>-t^'S,{x,\))x‘‘'-^‘-, (26) 

i=i 

s = 1, n. 

Proof. When s = 1 equality (25) follows from (23) for a; = 0, while (26) is obviously true. 
Suppose now that (25) and (26) have been proved for s = 1,..., iV — 1. Then 

N-l 

{(Uu,{x,p) - UUx.p)) - - 62,p«-'»)a^'‘'"'“S,(i,A)))i«-™ 

i=i 

N-2 

= ({UUx,p) - Ul^{x,p)) - Y,{hk,{p)p-^^ - A)))a;^i-^--ia;^--i-^- 

i=i 

-{hk,N-i{p)p~^^ - =J^kN{x,p), 

which gives (26) for s = N. We now write (23) as 

n 

T,m(x,p) = J2Mp)p~'‘' - 62iP"'-'“)p'‘'-™Sj(i, A) 

j=N 


+ 5^(t"y(pi)'‘'-"(4(i,A) - p(i,A))i«-™. 

J = 1 

Hence, using (17), we infer J7i:v(0,p) = {bkN{p)p~^^ — b\^p^^^~^^)cNo, which gives (25) for 
s = N. Lemma 4 is proved. 

Now we write (18) for z/ = 0 as 

J^u(x,p) = - r(f'(C/'i,(x,p)£;f («,p))(pi)"-‘-'‘“l4(i,p)rf« 

p ■'» j.l 



n 

U%{.x, p)Uj’*{t, p) = p^""'"W"^H^"”+^"p(a;, f. A), 
i=i 

it follows by way of (15) that 

n 

\'^U^o{x,p)Uj’*{t,p) < M3\{px)>^"->^^\, 0<t<a;<|p|"^ (28) 

i=i 
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Lemma 5. The following relations hold 

/ oo 

( Z l>%F;{pt)U°^‘{t,p)jVW,p)dt, (29) 

j=k+l 

J^ks{x,p) = (^'^^U^o{x,p)Uj'*{t,p)Ypt)'^-^->^-Vk{t,p)dt 

poo / ^ 

-Vi"-"- / I V b%^‘-i‘'(Mx,X)-Ct(x,\))F;(pt)uf’(t,p))vt(t,p)dt 

e=i j=k+i 

/ CXD ^ ^ 

(E (*.rt)i4(i,rt*), i<ipi-‘. (30) 

j=fc+l ^=s 

Proof. For s = 1, (29) and (30) follow from (27), in view of (22). Suppose now that (29) 
and (30) have been proved for s = 1,..., iV. Then, using (24), we calculate 

J^k,N+l{x,p) = {J^kN{x,p) - 7kN{^,p)SN{x,\)cflQ)x^^~^^+^ 

= -jfp-f - [ {^x%{x,p)uf-xp))(ptr-'-‘‘'v,{t.p)dt 

N—1 «oo / ^ 

- ^ b»,p“<->"(.%{x, Xj - C,{x, Xj)F;(pt)Up(f,pj)Vt(t,pj dt 

£=1 j=k+l 

/ oo ^ 

(5^(5^ b%p'"^~'"^x>^^->^^C^{x, A) - 65^p'^^-^^+i(7^(a;, A) 

j=k+l (=N 

-b%^«->‘«*‘(SK{x,\)-C,^(x,\))')F;{pt)Up{t,pm(t,p)dt) 

N n 

Y (S,(x, A) - C,(x, \))F](pt)Up{t, p)) V4(«, p) dt 

t=l “ ^ j=k+l 

/ oo ^ ^ 

( E ( E bxP^^-‘‘'x^’-^''^Xi(x,X))F;(pt)Up(t,p))vW,p)dt), 

j=k+l ?=v+l 

i.e. (30) is valid for s = N + 1. We now let a; —)■ 0 in (30) for s = iV + 1 , using (28), to obtain 
(29) for s = N + 1. This proves Lemma 3. 

It follows from (25) and (29) that 

1 poo XI 

h.(p)p-‘“ - bl =-^ (Y'’%F;(pt)Up(t,p)W,(t,p)dt. (31) 

P -^0 ^ j=k+l 

Using (31), (14), (19) and Lemma 3, we obtain 

hks{p)p~^^ - hi, = 0{J{p)) = 0(p"^), IpI ^ (50, p e Skp,, 


-E 


X‘ 


p-i—r-N+i 


i.e. (21) is valid. Theorem 2 is proved. 



3. Using constructed fundamental systems of solutions on each edge we can study spectral 
properties of systems on graphs and solve the inverse spectral problem. Let {ikj}k=T^j be the 
roots of the characteristic polynomial 

Uj fM-l 

^ ^ ('C ^rij,j ■ 1) — ■ 0. 

fi=0 k=0 

For dehniteness, we assume that ^kj — ^mj 7^ srij, s G Z, Re^ij < ... < Re ikj 7^ 0, rij — 3 
(other cases require minor modihcations). We set 6j := Uj — 1 — Re ~ (ij), and assume 

that the functions z/ = 0,/i —1, are absolutely continuous, and E L(0, Ij). 

Fix j = l,p. Let the numbers Ckjo, k = l,nj, be such that 

J_ J_ ^kjO ^b.et[(^^j ]k,u=l,nj^ 
k=l 

According to results of the previous section one can construct the fundamental systems of so¬ 
lutions {Skj{xj, of Eq. (1) on the edge ej such that the functions {xj, X), v = 

0,nj — 1, are entire in A, and for each hxed A, and Xj —)■ 0, Skj{xj,X) ~ Ckjox^^^ Consider 
the linear forms 

V 

UjviVj) = ^ljupyj^\lj), j = TTp, = 0,nj - 1, 

p=0 

where are complex numbers, := Ijuu 7 ^ 0. Denote (n) := (|n| +n)/2, i.e. (n) = n 

for n > 0 , and (n) = 0 for n < 0. Fix s = l,p, k = l,ns — 1. Let = {'4^skj}j=Xp be 
solutions of Eq. (1) on the graph T under the boundary conditions 

'ipsksixs, A) ~ Cksoxf^‘, Xs -E 0, (32) 

^skj{xj, X) = Xj ^0, j = T^p, j ^ s, (33) 

and the matching conditions at the vertex vq : 

Uiui'ipski) = Ujui'tpskj), j = ‘2,p, U = 0,k-1, Uj > u + 1, (34) 

p 

Uj^{^|Jskj) =0, u = k,ns-l. (35) 

j=l,nj >u 

Dehne additionally -ipsnsixs, A) := Sns{xs, A). Using the solutions {S^j{xj, A)} , one can write 

Uj 

'ipskjixj, A) = Mskjp{X)Sf,j{xj, A), J = l,p, k = 1, Us - 1, (36) 

p=i 

where the coefficients Mskjp{X) do not depend on xj. It follows from (36) and the boundary 
condition (32) for the Weyl-type solutions that 

ns 

^jjsksixs, X) = Sks{xs, X) + y] Mskp{X)S^s{xs,X), Mskp{X) ■.= Msksp{X). (37) 

^=k-\-l 

Denote 

M,(A) = [Mskp{X)]p^^^-Yjjp^, MskpiX) := 4^ for k > u. 

The matrix Ms{X) is called the Weyl-type matrix with respect to the boundary vertex Vg ■ The 
inverse problem is formulated as follows. Fix w = 2, p. 
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Inverse problem 1. Given {Ms(A)}, s = l,p\w , construct q on T. 

Fix s = l,p, k = l,ns — 1. Substituting (36) into boundary and matching conditions (32)- 
(35), we obtain a linear algebraic system with respect to Msfcj^(A). Solving this system one gets 
Mskjij,{^) = ^skW, where the functions and Asfc(A) are entire in A. In 

particular, 

Msk^J,{^) = Ask^{\)/Ask{^), k < p, 

where Ask^l{\) ■= Asksf^W- 

4. In this section we obtain a constructive procedure for the solution of Inverse problem 1 
and prove its uniqueness. First we consider auxiliary inverse problems of recovering the differen¬ 
tial operator on each each hxed edge. Fix s = l,p, and consider the following inverse problem 
on the edge e* . 

IP(s). Given the matrix Mg , construct the potential on the edge . 

Theorem 3. Fix s = l,p. The specification of the Weyl-type matrix Mg uniquely deter¬ 
mines the potential qg on the edge Cg . 

We omit the proof since it is similar to that in [1, Gh.2]. Moreover, using the method of 
spectral mappings, one can get a constructive procedure for the solution of the inverse problem 
IP{s) . It can be obtained by the same arguments as for n-th order differential operators on a 
hnite interval (see [1, Gh.2] for details). 

Fix j = l,p. Let ipjk{xj,X), k = l,nj, be solutions of equation (1) on the edge ej under 
the conditions 

= ^ku, ir =l,k, (Pkj{xj, A) = Xj 0. 

We introduce the matrix mj{X) = [mjku{X)]k^^=T^-, where mjkfiX) := (ppr^\lj, X). The matrix 
mfiX) is called the Weyl-type matrix with respect to the internal vertex uq and the edge Cj . 

IP[j]- Given the matrix mj , construct qj on the edge ej . 

This inverse problem is the classical one, since it is the inverse problem of recovering a higher- 
order differential equation on a hnite interval from its Weyl-type matrix. This inverse problem 
has been solved in [1], where the uniqueness theorem for this inverse problem is proved. Moreover, 
in [1] an algorithm for the solution of the inverse problem IP[j] is given, and necessary and 
sufficient conditions for the solvability of this inverse problem are provided. 

Fix j = l,p. Then for each hxed s = l,p\j, 


mjifiX) 


dy‘bj.A) 

fi’sijilj, A) 


z/ = 2, Uj, 


(38) 




det[figf,j{lj, A),..., A)] 


J fi=l,k 


, 2 < k < p < Uj. 


(39) 


Now we are going to obtain a constructive procedure for the solution of Inverse problem 
1. For this purpose it is convenient to divide diherential equations into m groups with equal 
orders. More precisely, let Ui > U 2 > ... > oJm > ^m+i = 1, 


Up. : Uj, 


j = l,m, 0 = po < Pi < ... < Pm = p. Take N such that pn = uo. 

Suppose that we already found the potentials qg , s = l,p \ pat , on the edges Cg , s = 
^jP\Pn. Then we calculate the functions Skj{xj,X), j = l,p\pAr; here k = l,Ui for 
j = Pi-1 + l,Pi. 


Fix s = l,pi (if iV > 1 ), and s = l,pi — 1 (if iV = 1 ). Our goal now is to construct the 
Weyl-type matrix mpj^{X). According to (38)-(39), in order to construct mp^(A) we have to 
calculate the functions 




( 40 ) 
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We will find the functions (40) by the following steps. 

1) Using (37) we construct the functions 

A), /c = 1 , Wjv - 1 , z/ = 0 , Wi - 1 , (41) 

by the formula 

A) = St\l„ A) + ^ hWWS'ZHl., A). (42) 

2) Consider a part of the matching conditions (34) on 4/^^ . More precisely, let ^ = N,m, k = 
u^+i,u^ — 1, I = ^,m, j = l,pi — 1. Then, in particular, (34) yields 

Upi^ui'fpskpi) = Ujv(^sfci), v = wz+i - 1, mm{k -l,ui-2). (43) 

Using (43) we can calculate the functions 

A), li = N, m, - 1. Z ^ ^ l.P!, ^ - l,min(fc - l.w; - 2). 

(44) 

3) It follows from (36) and the boundary conditions on that 

dgft.A)= 5 ; (45) 

/i.=max(LJ/—1) 

k = l,ui-l, / = l,m, j = pi_i + l,pi\s, u = 0,ui- 1 . 

We consider only a part of relations (45). More precisely, let ^ = N,m, k = u^+i,u^ — 1, I = 
l,m, j = pi_i + l,p/, j ^ Pn, j s, v = 0,mm{k 2). Then 

^ A) = 'ipillilj, A), z/ = 0, mm{k -l,ui-2). (46) 

/i=max(cj^—fc+1) 

For this choice of parameters, the right-hand side in (46) are known, since the functions (44) 
are known. Relations (46) form a linear algebraic system agkj with respect to the coefficients 
MskjpW- Solving the system askj we hnd the functions Mgkjpi^)- Substituting them into 
(45), we calculate the functions 

^£(^i>A), k = l,ujN-l, l = j =pi_i + l,pi\pN, Jy = 0,uji-1. (47) 

Note that for j = s these functions were found earlier. 

4) Let us now use the generalized Kirchhoff’s conditions (35) for . Since the functions 

(47) are known, one can construct by (35) the functions (40) for k = I^ojn — 1, = k^UN — 1- 

Thus, the functions (40) are known for /c = 1, ujq — 1, z/ = 0, Ujsi — 1. 

Since the functions (40) are known, we construct the Weyl-type matrix mpj^{\) via (38)-(39) 
for j = Pn- Thus, we have obtained the solution of Inverse problem 1 and proved its uniqueness, 
i.e. the following assertion holds. 

Theorem 2. The specification of the Weyl-type matrices Ms{X), s = l,p\pN , uniquely 
determines the potential q on T. The solution of Inverse problem 1 can he obtained by the 
following algorithm. 

Algorithm 1. Given the Weyl-type matrices Ms{\), s = l,p\pN ■ 

1) Find qs , s = l,p\pN , by solving the inverse problem IP{s) for each fixed s = l,p\pN ■ 

2) Calculate A), j = l,p \ ptv / here /c = 1, Wj, z/ = 0, — 1 for j = Pi_i -|- l,pi. 

3) Fix s = l,pi (if N > 1 ), and s = l,pi — 1 (if N = 1 ). All calculations below will he 

made for this fixed s. Construct the functions ( 4 I) via (4^)- 
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4) Calculate the functions (44) using (43)- 

5) Find the functions by solving the linear algebraic systems agkj ■ 

6 ) Construct the functions (40) using (35). 

1) Calculate the Weyl-type matrix mp^{\) via (38)-(39) for j = Pn ■ 

8 ) Construct the potential qp^ by solving the inverse problem IP[j] for j = Pn ■ 

We note that inverse spectral problems for Sturm-Liouville operators on spatial networks 
were studied in [2-7], and for higher order operators in [8]. 
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